Self field of ac current reveals voltage— current law in type-II superconductors 
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The distribution of ac magnetic fields in a thin superconducting strip is calculated when an ac 
current is applied to the sample. The edge barrier of the strip to flux penetration and exit is taken 
into account. The obtained formulas provide a basis to extract voltage-current characteristics in 
the bulk and at the edges of superconductors from the measured ac magnetic fields. Based on these 
results, we explain the spatial and temperature dependences of the profiles of the ac magnetic field 
measured in a Bi2Sr2CaCu208 strip [D.T. Fuchs et al.. Nature 391, 373 (1998)]. 

PACS numbers: 74.25.Qt, 74.25.Sv 



I. INTRODUCTION 



Recently a novel experimental approach for mapping 
the distribution of the transport current across a flat su- 
perconducting sample was developedi^ In this technique 
an array of microscopic Hall sensors is attached directly 
to the sample and the perpendicular component of the 
magnetic self-field generated by an ac transport current 
is measured at various locations across the sample. In- 
verting the Biot-Savart law, the distribution of the trans- 
port current across the sample can be obtained from the 
measured self field. It was found that in contrast to the 
common assumptions the transport-current distribution 
is highly nonuniform and varies significantly as a func- 
tion of temperature T, the applied magnetic field i?a, 
and the phase of the vortex matter. Fig. 1. Surprisingly, 
over a wide range of the H-T phase diagram the cur- 
rent flows predominantly at the edges of the sample due 
to significant surface barriers both in high- Tc and clean 
low- Tc superconductorsiiiSi^ In this situation the veloc- 
ity of the vortex lattice is determined by the rate of the 
vortex activation over the surface barrier rather than by 
the bulk vortex pinning. At lower temperatures bulk vor- 
tex dynamics takes over, resulting in redistribution of the 
transport current; see Fig. 1. So far, however, the experi- 
mental data thus obtained were not used for quantitative 
study of vortex dynamics due to the absence of a detailed 
theoretical description. 

In this paper we develop such a description. It provides 
the basis for the extraction of various voltage-current 
characteristics in the bulk and at the edges of supercon- 
ductors. In Sec. II we give general formulas for the ac 
magnetic fields and currents, while in Sees. III-V three 
simple models are considered which shed light on the 
data presented in Fig. 1. Combinations of these models 
(or of more complicated models developed on the basis 
of the general equations of Sec. II) enable one to analyze 
various experimental situations in different superconduc- 
tors. The obtained results show that the self-field method 
can be a useful tool for determination of the voltage- 



current characteristics of superconductors. This devel- 
opment supplements the usual transport measurements 
and thus opens novel possibilities for a quantitative in- 
vestigation of vortex dynamics. 



II. GENERAL FORMULAS FOR AC 
MAGNETIC FIELD AND CURRENT 

A. Edge barrier 

Consider an infinitely long thin superconducting strip 
of width 2w and thickness d ^ w, filling the space 
\x\ < w, \y\ < OO, \z\ < d/2. Let the strip be in a con- 
stant and uniform external magnetic field Ha directed 
along the z axis, and let an ac current la — lao smtot of 
frequency to and of magnitude lao be applied to the sam- 
ple. The current-voltage dependence for the supercon- 
ducting strip in the bulk is implied to have the standard 
form E = p(|j|)j, with resistivity p being a nonlinear 
function of the current density j : At small current densi- 
ties, \j\ < jc, the resistivity p (and thus the electric field 
E) is practically equal to zero; near the critical current 
density jc (i.e., at |j| w jc) the resistivity sharply in- 
creases, and it reaches a constant value pg independent 
of j in the flux- flow regime at \j\ :s> jc- 

In the narrow regions of the strip near its right and left 
edges X = where the barrier for flux penetration into 
the strip or exit from it occurs, we shall use the current- 
voltage laws: E = Rrlr, E = Rili where and /; are 
the currents that fiow in these edge regions while Rr and 
Ri are the appropriate resistances per unit length in y. 
In other words, we model the small edge regions of the 
strip as "wires" with the resistances Rr and Ri . Since the 
edge barrier is generally different for flux entrance into 
the sample and for flux exit from it, "4 these resistances 
are generally different, too. When la > 0, the vortices 
enter the strip at x = —w and exits from it at x = w, 
i.e., one has Ri — Ren and Rr = Rex, while at /a < we 
arrive at the opposite relations: Ri = Rex and Rr = Ren- 
in general. Rex and Ren are nonlinear functions of the 
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FIG. 1: Top; The temperature dependence of the in-phase 
part of the first harmonic of the ac self field measured by an 
array of 7 Hall probes at the surface of a Bi2Sr2CaCu208 
strip carrying an ac current of 4 mA amplitude in a con- 
stant applied dc field of 1000 G. Bottom: The temperature 
dependence of the second harmonic of the ac magnetic field 
measured by one of the sensors at 10 mA ac current at 1000 
G. Shown are also the characteristic temperatures Td, Tst, T^. 
AU the data from Ref. |3- The inset shows the geometry of 
the experiment and the positions of the sensors. 



currents flowing in the appropriate edge regions 

We may estimate the radius of the two edge wires from 
the following arguments. For a strip in a parallel mag- 
netic field when Ha is applied along the x axis and is 
larger than the penetration field, the z-size of the region 
where the surface currents flow is less than the London 
penetration depth It is generally believed that under 
condition A <C (i, which usually holds for single crystals 
of high-Tc superconductors, this is also true in the con- 
sidered case of a perpendicular magnetic field {Ha \\ z) 
if Ha is higher than the lower critical field Hd- In other 
words, one might expect that the width of the edge wires 
is less or of the order of A. However, simple considera- 
tions show that this is not the case since this assumption 
leads to a contradiction at least for not-too-large ratio 
Ha/Hci- Indeed, consider the situation when flux line 
pinning is negligible. If the edge currents flew only at 
w > \x\ > w — \, they would generate both H^ and H^ 
at least in the region w > \x\ > w — d. This means 
that the flux lines are curved at such x. But without 
pinning and currents the curved flux lines cannot be in 



equilibrium. Thus, we conclude that in a strip in a per- 
pendicular magnetic field the edge currents flow in the 
whole region w — < d and probably have a compli- 
cated distribution over x and z there. With this in mind 
we shall assume below that the width of the edge wires 
is of the order of d. Note that for Ha of the order of the 
lower critical field Hd this assumption agrees with the 
known result for the width of the geometrical barrier^ 



B. Simple case 

We now consider the magnetic field Hz{x, t) generated 
by the ac current la = laa sinwi. In particular, we calcu- 
late the first and the second harmonics of the ac magnetic 
field H^: 



Hi^Hx) = ^ f H,{x,t) sin ujtdt, 
^ Jo 

2 f'^ 

(2,) ^ / jj^ (2;^ ^) COS 2ujtdt, 
T Jo 



(1) 
(2) 



where T = 2tt/u! is the period of the oscillating ac cur- 
rent. (When the inductance of the strip is negligible, it 
follows from symmetry considerations that the first har- 
monic proportional to cos cut and the second harmonic 
proportional to sin2cL;i are equal to zero.) 

The total current la flowing in the strip at time t splits 
into the bulk current and the edge current /g, 

Ia=Ib+Ie, (3) 

where le consists of the left (/;) and right (7^) parts, 

Ie=Il+Ir- (4) 

If the frequency lo is not too large, the current distribu- 
tion over the cross-section of the sample is uniform, and 
hence the sheet current in the strip J{x) (i.e., the cur- 
rent density integrated over the thickness d of the strip) 
is independent of x. 



J{x) 



2w 



In this case we immediately obtain 

T - T T - T 

where Rb = p{J/d)/2wd and 

RlRr RenR: 

jXq — — 

Rl + Rr Re 

As to /; and Ir, one finds 

Rr 



Rp. 



Rb + Re 



Rk 



h = L 



Rl+Rr 



Ir — Ip, 



Rl 



Rl+Rr 



with 



Rl = Re 



Rr = Rex for la > 0, 



(5) 



(6) 



(7) 



(8) 



(9) 
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and 



General case 



Rl=Re 



Rr = Re 



for /a < 0. 



(10) 



As was mentioned above, R^n and R^x are gen- 
erally nonlinear functions of J = Ibl'^w and of the cur- 
rents len and /e:r, i.e., = Rbih), Ren = Renihn) and 
i?ea: = Rexihx), whcrC /e„ = // , /go; = Ir for /a > 0, and 

len = In lex = h for /a < 0. Therefore, formulas ((HJ - 
are in fact equations for /; and Ir, 

Rr{Ir)Rb{Ib) 



— 



Ib=Ia 



[Rlih] ^ Rr{Ir)]Rb{Ib) + Rr{Ir)Rl{Iiy 

Rl{Il)Rb{Ib) 

[Rl {h ) + Rr {Ir)]Rb{h) + Rr {Ir)Rl {h ) ' 
Rr{Ir)Rl{Il) 



'Rl{Il)+Rr{Ir)]Rb{Ib)+Rr{Ir)Rl{Il)' 



(11) 

(12) 
(13) 



On determining //, Ir, and lb from this set of equations, 
one can calculate {x, t) as sum of the fields of the wires 
and of the spatially uniform sheet current J, 



2tt \ 2w 



w — X 



w + X 



Ir 



W ~ X W + X 



.(14) 



This formula together with Eqs. Hll|) - (|13|) gives as a 
function of x and la, — IIz{x,Ia). Expressing the 
time t in Eqs. Q via la, sinwi = la/lao = la, we 
arrive at 



Hi'\x)-- 



- (\H,Xx,ia)+H^{x,-ia)] ^]-^2tn 

ttJq (1-4')'/' 



(15) 



(16) 



Note that II^{x,-Ia) ^ -II^{x,Ia) only if Ren ^ Rex- 

Hence, the second harmonic Hi {x) appears only when 
an edge barrier exists, and when this barrier leads to an 
asymmetry with /; ^ Ir- Using formulas (01, I©, p4|l . 
expressions H15|l . H16(l can be also rewritten in the form 



H^Pix)-- 



4:TTW 



w — x 



w + x 



w 



27r(w2-a;2) ^ ' 



+4'^^l (17) 



(18) 



where li^^ is the first harmonic of le, while Ale^' is the 
second harmonic of Ir — h. 



(2) . 







(1-/2)1/2 



Formulas (j5|l- H13|l have been obtained under the as- 
sumption that the resistances are much larger than ^qlo 
(quasistatic approximation) , 



Rl, Rr, Rb ^ Mow, 



(19) 



where Rb ~ p{J / d) /2wd. Under this assumption one can 
neglect the inductance of the sample, which is of the order 
of /^o (per unit length) Besides this, when Rb ^ A*ow, 
the two-dimensional penetration depthSiS A = pg/ /zgwd 
of the ac field in the Ohmic regime is considerably larger 
than the width 2w of the strip, 



2w pqlu 



(20) 



and hence, one may expect that the current distribution 
over the cross section of the sample is indeed uniform. 
Since the resistances Ri, Rr, Rb decrease with decreasing 
temperature T, the assumption (|19() is true for not too 
low temperatures. We now address the general case when 
conditions H19|l are not necessarily fulfiUedpi^ 

Let Ea be the voltage drop (per unit length) which 
just generates the current la- We now consider not only 
the edge regions but formally also the whole strip as a 
set of parallel wires. Then, we arrive at the following set 
of equations for the left and right edge wires and for the 
wires of width dx at points x {—w < x < w): 



Ea^RlIl+Lu,Il+LlrIr+ / dxJ{x,t)Lxi, (21) 

J —w 
nW 

Ea = RrIr-\-Ln,ir-\-Llril+ / dxj{x,t)Lj,r, (22) 



Ea^PxJ{x) + L^iIi^L^,rIr^ I dx J {x ,t)L^^' , (23) 

J —w 

Ia^Il+ Ir + h, (24) 

where px = p[J{x)/<I\/d is the sheet resistance, 

I dxJ{x,t), (25) 



the dot over currents means the time derivative, Lxx' = 
{p.o/2'K)\n{l/\x — x'\) is the mutual inductance per unit 
length of two parallel wires of length I located at points 
X and x' ^ Lxi ~ {po / 2tt) \n{l / \x + w\) and Lxr — 
(/io/27r) ln(Z/|a; — u;|) give the mutual inductances of the 
edge wires and that located at x, Lir = {po/2tt) ln.{l/2w) 
is the mutual inductance of the two edge wires, and 
Lw = (Mo/27r)[ln(Z/r) -1-1/4] is the self inductance of the 
edge wires. Here we have assumed that both edge wires 
have the characteristic radius r ^ d/2. The magnetic 
field Hzix, t) is given by the Biot-Savart law 



H,{x,t) 



Ir 



h 



1 



27r \w ~ X w + x 



"J{x',t) dx'' 



(26) 
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Note that under conditions (|19|l . one may keep only the 
first terms on the right hand sides of Eqs. (|21|I - H23|) (omit- 
ting all time derivatives) and replace formula H25|l by 
Eq. JSj) (since the skin effect is absent). After this sim- 
plification, equations (|21|) - H24|) and H2()l) are equivalent 
to formulas © - (HHl. 

In some region above temperature Td shown in Fig. 1, 
one has Rt ^ Re- Thus, a situation is possible in which 
Rb ^ /io'jJ > Re, and the skin effect is negligible, but one 
cannot omit completely the terms with the inductances 
in Eq. (|21|I - H23|) . In this context, it is useful to consider 
the case when the only assumption is the uniform distri- 
bution of the sheet current over x, i.e., the fulfilment of 
Eq. (0) . Under this assumption Eq. H21() - H24() reduce to a 
form which describes three parallel connected conductors 
with inductive coupling: 

Ea=RlIl+L.ujIl+LirIr+LwbIb, (27) 

Ea=RrIr + LwIr+LlrIl+Lwbh, (28) 

Ea=RbIb+LwbIl+LwbIr+LbIb, (29) 

la^Il+Ir+h, (30) 

where Rb = p{j)/2wd with j = Ib/2wd, Lb = 
(/Lto/27r)[ln(Z/2w) -t- 3/2] is the inductance of the strip 
for uniform current distribution in it, and L^jb — 
(/Lto/27i")[ln(//2w) -I- 1] is the mutual inductance of the 
strip and one of the edge wires. In this case Eq. 
transforms into expression H14(l . and formulas H17|l . IjlSI) 
are valid. 

Equations H21|) - (|24|l generalize the known equation for 
a strip with time-dependent current p^S which follows 
from Eq. H23|) if one omits the edge- wire terms with Ii and 
jr. Thus, the numerical procedure elaborated in Refs. 0, 
13 can be applied to solve Eqs. H21|) - (|24|) in the general 
case. However, in this paper we shall confine our anal- 
ysis to simple models which are sufficient to understand 
the main features of the experimental results obtained in 
Refs.EIll 



III. OHMIC MODEL 

We begin our analysis with the simplified symmetric 
Ohmic model in which Ren — Rex = 2_Re, and these re- 
sistances are independent of the currents len — lex but 
depend on temperature T. We shall use a temperature 
dependence that satisfies the following requirement: At 
T = Tc when the edge barrier is absent, the resistances 
are proportional to Rb with a geometrical factor that is 
the ratio of the cross-section area of the strip to the cross- 
section area of one of the edge wires. At lower temper- 
atures we then take the conductance of the edge wires, 
1/ Rem ^/Rex, as a sum of this bulk contribution and the 
conductance caused by the edge barrier: 



1 



1 



1 



ReniT) RexiT) i?o(r) 2R0 



1 f Ue{t) 

exp— 1 , (31) 



where Ro{T) — {w/r)Rb{T), i?g is some constant, t = 
T /Tc, Tc is the temperature of the superconducting tran- 
sition, Ue{t) is the dimensionless magnitude of the edge 
barrier (in units of Tc), and we shall imply a linear tem- 
perature dependence of this barrier, Ue{T) = ?7e(0)(l— t), 
with some constant ?7e(0). The function Ro{T) just gives 
the resistance of one of the edge wires with the edge 
barrier neglected (r ^ d/2 is the characteristic radius 
of the edge wire), while the second term in Eq. (|31|) 
is the edge contribution described by an Arrhenius law. 
We have introduced the constant —1 in the brackets to 
take into account that this part of the conductance van- 
ishes at Tc. Besides this, in this Ohmic model we use 
jc{T) — 0, and thus Rb is independent of J at all tem- 
peratures. For definiteness we take the temperature de- 
pendence of Rb in the form of an Arrhenius law, too, 
R^{T) = pff/2wd = Rb{Tc) exp[-Ub{0){l-t)/t], but with 
the constant Ub{Q) smaller than Ue{0). The requirement 
Ub{0) < Ue{0) emphasizes that with decreasing temper- 
ature the main contribution to the conductance of the 
edge wires results from the edge barrier. In the discus- 
sion below, we shall imply parameters that correspond 
to the experimental valuesi^ lao = 4 mA, w — 75 pm, 
d = 10 iim, 2w/d = 15 (f* w/r), I = 1.5 mm, Tc = 88 
K, a;/27r = 73 Hz, Rb{Tc) « 6.7 n/cm, yielding the large 
value A/w = 2.32 • 10® for the two dimensional penetra- 
tion depth at T = Tc- For definiteness, in the construc- 
tion of figures we shall also take R^ = Ro{Tc). As to 
Uc{0) and Ub{Q), we chose them so that one can repro- 
duce the characteristic temperatures T^ and T^b marked 
in Fig. 1. 

Within this simplified model, one has Hi^\x) = 

(2) 

Hz{x,t)/ sin ujt and Hz (x) = 0, where under conditions 
H19|l (at sufficiently high temperatures) Hz{x,t) is given 
by Eq. (|14() . while Ie{= 2Ii = 2/,.) and lb are described 
by Eqs. ©• Thus, we obtain 



HPix)^^(. 



Re 



A-KW \Re + Rb 



In 



w+x 



Rb 



2wx 



Re + Rb 



(32) 



With decreasing temperature when the inequality Re < 
jjLQOj becomes valid, it is necessary to take into account 
the inductances of the strip and of the edge wires. In this 
case Eqs. ()27 |l - (|30|) reduce to a set of linear equations (the 
time derivatives are replaced by the factor icj), and after 
simple calculations we find 

2Ii = 2Ir=Ic=Ia^^, Ib = lazr^, (33) 

Zb + Ze Zb + Ze 

where Zb and Zg are the "effective impedances" of the 
strip and of the two edge wires, 

Zb^ Rb + i(^{Lb- L^b) = Rb + ii^^ , (34) 

Ze^ Rc + iuj(^ +^ -Lyjb 



fiQ f 2w 7 

Re + In - 

47r \ r 4 



(35) 
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Note that even though the various inductances calculated 
per unit length of the strip depend logarithmically on its 
length I, the resulting effective impedances Zb, (per 
unit length) are independent of I. Inserting expressions 
(|33|l in formula (|14|l . one obtains Hi^\x). In other words, 
down to the temperature Td defined in the Ohmic model 
by Rb{Td) ^ /^o^, the first harmonic is described by for- 
mula H32I) in which i?f, and Re are simply replaced by 
Zfc and Zg. The real part of this expression gives the 
harmonic which is in-phase with /„, i.e., proportional to 
sinwt like lait), while the imaginary part corresponds to 
the out-of-phase signal proportional to cos ujt. 

In Fig. 2 we show the temperature dependence of 
KeHz^^ and hnHz^^ for discrete values of x similar to the 
locations of Hall sensors in the experiment We first 
discuss the temperature dependence of the in-phase sig- 
nal TleHz^\ At temperatures near Tc, the bulk current lb 
exceeds le since Rb{Tc) < Re{Tc), and one has an Hz[x) 
that is characteristic of the normal state of the supercon- 
ductor [the first term in Eq. (|32|l dominates], see Fig. 3. 
But the ratio of the currents I^, and lb changes with tem- 
perature. The role of the edge current increases with 
decreasing T since C/f,(0) < C/e(0), and at a temperature 
T = T^b « 69 K defined by Re{Tsb) = Rb{Tsb) one has 
lb = le, and the contributions of and /g to Hz practi- 
cally compensate each other for x inside the strip. With 
further decrease of T, the x-dependence of Hz is mainly 
determined by the edge current [i.e., by the second term 
in Eq. 1(221)] ■ When the inductances begin to play a role 
at temperature T^; w 46 K defined by Re{Tx) ^ Hqw [or 
more exactly, by Ileze(Tx) ~ lmZe{Tx)], the impedance 
Zf, saturates while Zb still decreases with decreasing T. As 
the temperature ~ 38 K corresponding to Rb{Td) ~ 
Hquj is approached, the contribution of to Hz begins to 
increase again, and the ratio of the currents, Ib/Ie tends 
to the constant ln(2u;/r) — 7/4, which is determined by 
the ratio of the imaginary parts of Ze and Zb- Of course, 
in reality at such temperatures the skin effect plays an 
important role, and ReTji"'^' tends to zero rather than 
to a profile determined by the above-mentioned ratio of 
the currents, see Fig. 2. As to the out-of-phase part of 
the first harmonic Im7ji^\ it appears only near the tem- 
perature Tx at which the inductances begin to play a 
role, and it differs from zero when the real and imagi- 
nary parts of the effective impedances are comparable. 
At still lower temperatures when Re and Rb become neg- 
ligible as compared to the imaginary parts of Ze and Zb, 
the out-of-phase signal vanishes again. 

We now present the results for the asymmetric case 
when Rex ^ Ren- In this situation model (|31|l is gener- 
alized as follows: 



1 



Rex {T) 
1 



Ra{T) 
1 



1 

w 

1 

W 



exp' 



cxp- 



Ue{tl 
t 

Ue{t) 



1 



- 1 



(36) 



Ren{T) Ro{T) ■ i?u„ V"""^ t 
where i?gj,(T), i?g„(r) are some constants, and Ro{T) is 
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FIG. 2; The temperature dependence of the first liarmonic 
of the ac self field which is in-phase with the applied cur- 
rent, KeHi^\ (top) and out-of-phase with it, lmHi^\ (bot- 
tom) at positions —x/w = 1.2, 0.8, 0.6, 0.4, 0.2, on the 
surface of a long thin strip of width 2w. The applied ac 
current is la ~ lao sin cot. Results within the symmetric 
Ohmic model with Ren = Re:^ for Ue{0) = 21, Ut{0) = 12, 
Rt,{Tc) = 6.7 n/cm, w/r = 15, R° = 15Rb{Tc), uj/2tv = 73 
Hz, Tc — 88 K, w = 75 fim, lao = 4 mA. The dashed lines 
are the analytical results, Eq. 138L while the solid lines are 
numerical results from Eqs. 1211 - 126(1 allowing for the skin ef- 
fect. The solid and dashed lines coincide at T > 42 K, see 
the enlarged plot (middle). Only numerical results are shown 
for ImHi^^ (bottom) since the analytical results, Eq. 13811 . ne- 
glect the skin effect and thus are applicable only at T > 42 K. 
Note the different scales of the T axis for the different plots. 
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FIG. 3: The in-phase components of the magnetic field and 
of the sheet-current profiles, Hi^\x) (solid lines) and J{x) 
(dashed lines), for the symmetric Ohmic model {rj = 0) at 
temperatures T = 85, 70, 50, 38, and 30 K, accounting for 
the skin effect which occurs at T < 42 K, cf. Fig. 2. The 
parameters of the model are as in Fig. 2. The circles mark 
the positions of the Hall probes used in Fig. 2. Both Hi^^ 
and J are in units of lao/iv. The total current in the depicted 
half strip is thus 1/2, equal to the area under the dashed lines 
plus the (not shown) edge current at a; = to. 



the same as in Eq. H31|l . It is convenient to introduce the 
parameter 



'0 
ex 



R 





en 



(37) 



which characterizes the asymmetry of the edge barrier. 
Then, one has 

R",„^2R"J{l + r^), = 2i?o/(l-ry), 

where i?^ = R°„R°^/{R°„ + Rex)- We emphasize that 
the constant rj defines the ratio {Rex — Ren) /{Rex + Ren) 
at sufficiently low temperatures when the conductance 



1/Ro is negligible in Eqs. while at T ^ Tc the 

resistances Ren and Rex tend to Rq, the resistance of 
the edge wires without barrier. In other words, formulas 
(|36|l describe both the asymmetry of the edge wires at low 
temperatures, and the disappearance of the asymmetry 
when the edge barrier vanishes. Note also that Re = 
Rex Ren/ {Rex + Ren) is uot equal to i?^ but approaches 
it with decreasing temperature. 

In this asymmetric case the first harmonic becomes: 



where 



laO 



AWW \Zeri + Zh 



■In 



w~x 



w + x 



Zeri + Zb W' —X' 



{Rex - Ren)"^ ip 



A{Rex + Ren) 1 + ' 

Ze, Zb are given by formulas (|34(l . 1)35(1 . and 



p 



Re 



Re 



tt{R. 



IX H~ Ren) \ ^ 



(39) 



(40) 



This p characterizes the relative contribution of the in- 
ductance of the edge wires to their effective impedance Ze- 
Analysis of formulas H38|) - (|40|l shows that the parameter 
rj practically has no effect on i/i^'' within the considered 
Ohmic model, and hence the first harmonic in the asym- 
metric case, in fact, coincides with the harmonic in the 
symmetric case of 77 = 0. However, for the case rj ^ the 
second harmonic appears, 
is 



This hP{x), 



Eq. 



21, 



aO 



W 



{Rex Ren) 



Zb 



1 



i^Rex ~^Ren) Zb + Ze,, 1 + ip 



(41) 



The real part of this expression gives the harmonic pro- 
portional to cos2ti;t, and its imaginary part to sin2a;i. 
In Fig. 4 we show the temperature dependence of \Hz \ 
at X = 0.7w. At temperatures near introduced by 
ReZe{Tx) ^ linZe{Tx), one has p{Tx) ~ 1, the amplitude 
of the second harmonic begins to decrease with decreas- 
ing T and becomes small even at temperatures above 
Td if Rb{Tx) 3> /iow. Of course, at temperatures when 
Rb{T) ~ ^J■oUJ, the second harmonic Hz vanishes in any 
case due to the skin effect. 

We emphasize that the formulas of this section are 
valid not only for the model dependences H31|) or (|36l) 
used here but also for any functions Rb{T), Ren{T), 
Rex{T). Thus, one can obtain some information on these 
functions fitting the theoretical dependences of Hi^\T) 

and Hz (T) to appropriate experimental data. Indepen- 
dent information on the functions Rb{T), Re{T) for each 
specific sample can be also obtained from usual transport 
measurements.^^ Beside this, some data on this subject 
can be extracted from frequency dependences of H^\t) 
and H)r'{T) since these harmonics are functions of the 
combinations {Rb/noto), {Ren/p^Q^), {Rex/p^i)^)- In the 
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FIG. 4: The temperature dependence of the magnitude of 
the second harmonic of the ac self field, |-ffi'^''|, at position 
\x/w\ = 0.7 on the surface of a long thin strip of width 2tj; 
with applied ac current la ~ laosinut. Results within the 
Ohmic model, Eq. 14111 . for the same parameters as in Fig. 2 
but with asymmetry = 0.1 and amplitude lao = 10 mA 
(line at u>/2tt = 73 Hz). The other lines are for different 
frequencies. 

framework of our model H3fci|) , if one changes the frequency 
u from wi to UJ2-, the temperature dependence of the first 
harmonic changes only near T^, while the temperature 
dependence of the second harmonic changes in the vicin- 
ity of Tx (the appropriate parts of these dependences 
shift to higher temperatures with increasing lo), see, e.g., 
Fig. 4. The shifts of the temperatures and Tx are 

and they give information on Ub{0) and C/e(0). 

We now briefly discuss the applicability of the consid- 
ered temperature dependences of Re and Rb to BSCCO 
samples. The experimental data on resistance of BSCCO 
strips obtained from usual transport measurements are 
well approximated by an Arrhenius Iaw^iii24i^ii4 In prin- 
ciple, parameters Ue (0) and R^ can be found from these 
data, while Ub (0) can be estimated in similar experiments 
with wide samples.^^ In particular, in Ref.llllit was found 
that Ue{0) !=a 12 for the same crystal as in Fig. 2, and it 
was estimated that Ub{0) < 8. Both these values are no- 
ticeably less than those used in the construction of Figs. 2 
and 4. The latter values were chosen such that one can 
reproduce the experimental Td and Tsb- In Figs. 2 and 4 
we have also used the specific value of R^, R^ = Ro{Tc). 
Other choices of the parameter i?^ enable one, in princi- 
ple, to use a smaller value of Ue{0) — f76(0), but do not 
provide the coincidence of both C/e(0) and J7b(0) with 
their experimental values [in the Ohmic model the value 



of /7f,(0) w 12 is fixed by the condition Rb{Td) « fJ.oUj]i^ 
Such disagreement seems to signal that the Ohmic model 
does not work in the whole interval from to Tc for the 
BSCCO crystals. 

Finally, we note two characteristic features of the 
Ohmic model studied here. First, if the magnetic field 
is measured in units of lao/w, the first and the sec- 
ond harmonics are independent of the amplitude lao 
of the applied current. In particular, the characteris- 
tic temperatures T^, Tbs, Tx do not shift with a change 
of lao- This scaling property may be used in experi- 
ments to find the temperature regions where the Ohmic 
model is applicable. Second, the resistance of the strip, 
R = RbRe/ (Rb + Re) is also independent of lao and can 
be used for the same aim as well. 

IV. MODEL OF NONLINEAR R,{Ie) 

The data of Ref. '3 obtained for different currents lao 
show that at least for some superconductors Re cannot 
be considered as a quantity that is independent of le for 
all temperatures. To get an insight into this problem, we 
shall now analyze a simple model of nonlinear dependence 
Re{Ie) and show that this nonlinear! ty can provide an 
alternative description of the experimental temperature 
dependences of and H^'^ near Tsb- In principle, this 
and a similar nonlinear model for Rb{Ib) near Td enable 
one to describe the experimental data on and h'^'^ 
with realistic values of [4(0) and C/h(0). 

We shall consider the following simple model depen- 
dences Ren{Ien) and Rex{Iex) for the flux entrance into 
the strip and the flux exit from it: When I en (or lex) 
exceeds the critical current I en {lex) fo'' Anx entrance 
(exit), the resistivity of the edge wire is the same as in 
the bulk of the strip, and hence the resistance of the edge 
wire Ren{Ien) [Rex{Iex)] coiucidcs with Rq = {w/r)Rb. 
At len = I en i^ex = lex) tbc resistance sharply drops 
down to a value Ren (Rex) which is independent of the 

current len [I ex) for len < len i^en < Ien)i i'^-, One haS 
Ren(Ien) = Ro for ^en > I em 

Ren{Ien) — Ren for len < I em (43) 

and 

Rex(Iex) — Ro for lex > lexJ 

Rex (lex) = Rex for lex < I^x' (44) 

In other words, we replace the complicated nonlinear de- 
pendence of the edge barrier on the current flowing in 
the edge region^ by a sharp jump of the appropriate re- 
sistance. Although with the equations of Sec. II, the 
currents and magnetic field can be found for any realistic 
dependences Ren{Ien) and Rex{Iex), within this simple 
model the calculations can be carried out analytically. 
This enables one to get insight into situations with var- 
ious relations between the edge critical currents /g„, /g„ 
and the resistances Rem Rex- 
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FIG. 5: Top: Temperature dependence of the first harmonic 
Hi^^ of the ac self field at positions x/w = —1-2, —0.8, 
—0.6, —0.4, —0.2, on the surface of a long thin strip of 
width 2w with applied ac current la = lao sin uit. Results 
within the nonlinear iie-model for the same parameters as 
in Fig. 2 (symmetric case, rj — 0) but with f/e(0) = 12, 
Ut{0) = 2, /e"„(0)/7,o = IUO)/Iao = 4. Assumed was 
IUT)/IeAO) = IUT)/IeAO) = (1 - tf where t = T/T,. 
Bottom: Temperature dependence of the first harmonic li^^ 
of the edge current. The parameters are the same as in the 
top panel but /en(0)//aO = 4 (symmetric case), 5, 10, 100. If 
ri differs from zero, the curves remain practically unchanged 
for 77 < 0.3. For comparison, the dashed line shows lji^\T) 
for the symmetric Ohmic model with parameters of Fig. 2. 



Let us specify the details of this model. With de- 
creasing temperature the currents /e„ and I^^ increase, 
and for definiteness we take /g„ = /g„(0)(l — t)"^, I^^ = 
/g2,(0)(l — t)^ where t = T/Tc (such temperature depen- 
dence of Jg„-|-/g^ was obtained in Ref.0). We also imply 
that near the temperature Tgb these critical currents of 
the edge wires become larger than the amplitude of the 
applied current lao, and thus the drop of edge- wire re- 
sistances occurs at sufficiently high temperatures where 
conditions l(T^ arc true, and so the inductances and the 
skin effect are ignored below. As to Rem Rex and Rb, we 
choose the same expressions for them as in the Ohmic 
model, but now we use Ue{0) ~ 12 and Ub{0) ~ 2 (in 
units of Tc), which are essentially less than those used in 
Sec. III. 

To find the first and second harmonics of Hz, one 
should solve Eqs. Hll|l - H13|) and find the currents Ir, 



and as functions of la- From these currents, one cal- 
culates Hi^\ h'^P using formulas lfT7|l - lfTH|) . Within the 
above simple model, these calculations can be done ana- 
lytically. For example, in the symmetric case when rj — Q 
and Jg„ = = /e/2, one has formulas © for lb and 
= 2Ii = 2Ir at la < h = It ■ {Rb + Re)/Rb whcrc 
Rb and R^ — RenRexl {Ren + Rex) are independent of lb 
and /g. In the interval h < h < h = le ■ {'^ + w/2r) 
we find le = le, h — la ~ le, whilc at I2 < la equations 
(HU-lO yield 



Rb 



Rb + 0.5i?o 
0.5i?o 



/ w 



Rb 



= 1 - 
0.5i?o " 2r 



1 + ^)" 
2rJ 



(45) 



Using these formulas, wc find the following expression for 
the first harmonic of /„: 



2/aO Rb 



COS 02) + 



Rb + Re 
2/aO 



sin 20] 



4/: 



2 y TT 
sin 202 



(cos 01 



(46) 



where 0i = min(arcsin[/i//ao], 7r/2) and 02 = 
min(arcsin[J2//ao], 7r/2). Inserting Eq. (|46|l into formula 
H17|l . we obtain Hi^\x). The second harmonic of 
in the symmetric case is equal to zero. In Appendix A 
we present the first and the second harmonics of in 
the asymmetric case when rj ^ 0, /g„ 7^ 1^^. Note that 
within this nonlinear Re model the characteristic tem- 
perature Tsb is still defined by the condition /i^^* = ^ao/2 
as in the Ohmic model. 

In Fig. 5 we show the first harmonics Tji^-* of Hz{t) 
calculated within this nonlinear Re model. The temper- 
ature dependence of the first harmonic is qualitatively 
similar to that of Fig. 2, and with increasing asymmetry 
of the edge barrier it does not change essentially. On 
the other hand, the second harmonic is highly sen- 
sitive to the asymmetry, see Fig. 6. Figure 6 shows both 
the change of with /e«(0)/4^:,(0) at fixed 77 and its 
change with 77 at fixed /g„(0)//g^(0). Interestingly, the 

(2) . 

temperature dependence of Hz is nonmonotonic and 

looks qualitatively similar to that of Fig. 4. But now the 

f 21 

temperature where Hz (T) sharply decreases is due to 
the temperature Ti of Appendix A at which some char- 
acteristic current of the edge barrier becomes larger than 
lao- This temperature Ti is independent of w within our 
nonlinear Re model. At temperatures T <Ti during the 
whole period of the ac field, the currents flowing in the 
edge wires are less than the appropriate critical currents 

/g^, the strip is in the Ohmic regime, and hP{T) 
(2) 

and Hz (T) obtained from the formulas of Appendix A 
coincide with expressions H38|l and H41|) if one neglects the 
inductive terms proportional to fiQUj in these expressions. 
This coincidence permits one to combine this nonlinear 
model with the Ohmic model of Sec. Ill and thus to ob- 
tain continuous curves in the entire temperature region. 
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FIG. 6: The temperature dependence of the second harmonic 
of the ac field defined by Eq. (|5J at \x/'w\ — 0.7 within the 
nonlinear Re model, Eqs. p8ll . IIA2^ . The parameters are the 
same as in Fig. 2 but with Ue{0) ~ 12, Ub{0) = 2, and besides 
this lex (0) /lao ~ 4. The top part shows the second harmonic 
at fixed 77 = for diff'erent I^„{0)/Iao = 4.5, 5, 7, 14, 40, while 
the bottom part gives this harmonic at fixed /en(0)//aO = 6 
but for different asymmetry rj = 0.025, 0.035, 0.05, 0.1, 0.2. 
Note that H^^^ (T) for 77 / does not reach zero after its drop 
in the vicinity of 60 K. 

Interestingly, when T is larger than the temperature T4 
of Appendix A, the maximum currents in the edge wires 
exceed /g„ and I^^ , and the strip is in the flux- flow regime 
during some part of the ac period. In Fig. 6 this temper- 
ature T4 corresponds to the point where Hz (T) crosses 
zero (i.e., r4 « 78 K in Fig. 6b). 



p depends on the current density j. We shall analyze the 
following model current-voltage dependence: 

E^pU)j^Pff .^l^'^^ly (47) 

where 

e(j) = 2e-^o/^sinh(^^), (48) 

and jc and the effective depth of flux-pinning well Uq 
are some decreasing functions of T. Let us consider the 
case Uq/T :s> 1. Then, at j > jc one has e(j) » j/jc, 
and thus p{j) — ps, while at j < jc the quantity e(j) is 
small, e(j) < j/jc, and we arrive at p{j)j = psjc^ij) < 
Thus, equations (|T7jl . lead to a sharp jump of 
resistivity at j = jc that is characteristic for the critical 
state model. In analytical calculations we shall simplify 
the model dependence further, putting 

P(J) = Ps for j > jc, 

p(j) = for J < Jc. (49) 

As to the flux-flow resistivity pff , for deflniteness we shall 
imply the same temperature dependence as in Sees. Ill 
and IV, pff(T) = ps(Tc)exp[-Ub{0){l - t)/t], where 
Ub(0) is a dimensionless constant and t = T/Tc- Since 
in this section we are mainly interested in understand- 
ing the effect of nonlinear Rb{Ib) on the first harmonic 
of the ac field near Td, we assume that near Td one 
has i?fr = ps/2wd 3> p-ooj. Re, and the resistances 
Ren = Rex — 2i?e are independent of len — lex (i-e., 
we treat only the symmetric situation and imply that 
the skin effect occurs sufficiently below Td). The tem- 
perature dependence of the critical current Ic = 2wdjc 
is assumed to be a monotonically decreasing function. 
From our analysis it will be clear that this dependence 
plays the most important role in the vicinity of the tem- 
perature Tj at which Ic{T) reaches lao, and so we may 
use the following representation of Ic{T): 



V. MODEL OF NONLINEAR Rb{Ib) 

As was mentioned above, in the Ohmic model below a 
temperature Td defined by the condition RbiTd) ^ IJi-quj, 
the skin effect occurs, which means the perpendicular ac 
magnetic field Hz is expelled from the strip. However, 
there is another reason for this field to be expelled from 
the superconductor. Indeed, if the ac current fiowing in 
the bulk is smaller than the critical current Ic = 2wdjc of 
the strip, the ac field cannot penetrate completely into 
the sample. Thus, a finite Ic and its increase with de- 
creasing temperature can, in principle, explain the exis- 
tence of the temperature Td in Fig. 1, below which Hz is 
zero in the strip and at its surface. To describe this, we 
now consider a nonlinear Rb model in which the critical 
current density jc is not equal to zero, and the resistivity 



IciT) = laO 



Tc~T 



(50) 



where n is some exponent. Note that in this represen- 
tation a change of /ao requires a change of Tj such that 
IciT) remains unchanged. 

Let us first show that at lao ^ Ic — 2wdjc our non- 
linear Rb model reduces to the Ohmic model considered 
in Sec. III. This situation always occurs at sufficiently 
high temperatures when Ic is small. The nonlinearity of 
Rb{Ib) gives only small corrections to the Ohmic results, 
and these corrections can be calculated analytically. In- 
deed, at times t when /fc(t) > Ic, the Ohmic model is 
applicable, and we have 

h = -^60 sui{ujt + Lpb), le = leo sin{ujt + ipe) (51) 
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FIG. 7: The temporal dependences of the apphed current 
la = I ao sin u>t, the bulk current and the edge current 
Je obtained numerically from Eqs. 12H - l25j with the realis- 
tic current- voltage dependence 14711 . 1481 1 at T = 36.5 K and 
Uo/T = 30. The other parameters are the same as in Fig. 2 
but Ub{0) = 6, Ue{0) = 15, and ho = 3.577c. At this temper- 
ature formula 15611 gives h/Ic ~ 3700. With increasing lao 
the width of time intervals where It changes from Ic to —7c 
or vice versa decreases. 
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FIG. 8: The temperature dependence of the phase shift (pb 
of the bulk current Ib{t) as compared to the applied current 
Ia{t) = laosinuit, Eq. 1531 . The parameters are the same as 
in Fig. 2 but Ub{0) = 6, Ue{0) = 15. The temperature « 39 
K is defined by ReZe(7a;) = lmz^{Tx), i.e., by ipb ~ 7r/4. 



where Iw and /eo are the amplitudes of the currents lb 
and /e, 



^60 — laf) 



|^e| J |^e| 



Rtf 



(52) 



the effective impedance of the wires Zg is given by 
Eq. H35|) . and ift, fe define the phase shifts of and 
Je with respect to la = laO suiut, 



tan ipb 



ImZe 

Reze 



tan Lfe 



ImZe 



< 1. 



(53) 



Formulas H51|l - H53|) are valid for those times t when 
\Ib{t)\ defined by these formulas exceeds Ic. When the 
absolute value of Ib{t) reaches Ic, the bulk current 
remains equal to ±/c, while /g =Ia — Ib, see Fig.|7|''^i It is 
important that is uniformly distributed over the width 
of the strip, J{x) = ±Ic/2w = ±jcd, so long as — 
Ic- This is due to the sharpness of the current-voltage 
law H49() . which prevents the skin effect. A nonuniform 
distribution of the sheet current J over x occurs only 
in a relatively narrow time interval (as compared to the 
period T = 2n/uj) when lb changes from Ic to —Ic or vice 
versa. The relative width of this interval is determined by 
the small ratio Ic/Iao ^ 1, and in the first approximation 
we can neglect this time interval in calculating the first 
harmonic of Ib{t), 



2 f'^ 

/^'^ = - / hit) sin Lotdt, 
^ Jo 



(54) 



which is in-phase with the applied current. Then, we find 



21, 



60 



■ COS ipb ■ 



sin 2o92 2/9 , X , s 
-^ + -i(l-cos(p2) ,(55) 



where ip2 = min(arcsin[/2//ao], 7r/2), h is the total cur- 
rent at which the amplitude of the bulk current, ho, be- 
comes equal to h, 



J -7 

72 = 7c f-; 



Rs 



C I ~ 



(56) 



and Ibo, fb are determined by formulas If 
laO 3> I2, one has (^2 ^ 1, and formula (|55|l indeed re- 
duces to the result for the Ohmic model, I^^^ ~ ho cos ipb- 
On the other hand, at I2 > ho 3> h we have (^2 = '^/2, 
and this formula gives 

4^'' ~ -h cos ipb. (57) 

TT 

Thus, in this case the temperature dependence of /^^^ is 
straightforwardly expressed via the temperature depen- 
dence of the critical current and of the phase shift ipb- 
The dependence <fb{T) is shown in Fig. 8. Taking into 

account that li^"^ = ho — and using formula H17|l. 
which is always valid in the absence of the skin effect, 
we obtain the first harmonic of the magnetic field un- 
der condition lao ^ h, i-c, when the calculated H^z'\t) 
does not deviate considerably from the appropriate re- 
sult of the Ohmic model, see Fig. 9. Here H^^^ is the 
first harmonic defined by Eq. i.e., its in-phase part. 

In the nonlinear Rbih) model considered here the tem- 
perature Td is of the order of Tj, and it does not co- 
incide with the temperature at which Rb{T) ^ fJ-o^^, 
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FIG. 9: The temperature dependences of the first har- 
monic Hi^^ defined by Eq. Q within the nonlinear Rb model. 
The parameters are the same as in Fig. 2 but Ub{0) = 6, 
Ue{0) = 15. The critical current Ic is modelled by the power 
law 7c//ao = [{Tc - T)/{T, - Tj)]'^° with Tj = 37 K (top) 
and Tj = 40 K (bottom); Tc = 88 K. This shift of Tj from 
37 K to 40 K corresponds to the increase of Ic{T) by ap- 
proximately 6.16 times [or equivalently, to the decrease of lao 
by 6.16 times at fixed Ic{T)] and leads to the change of Td 
from Td « 38.5 K to Td « 42 K. The sohd fines show the 
numerical result obtained from Eqs. I2H - II25L 114 7L 11481 with 
Uo/T = 30, while the solid lines with dots are the analytical 
result, Eq. 15511 . The horizontal dashed lines mark //i^' in the 
Ohmic model when IciT) = 0. The other dashed lines show 
the analytical approximation obtained from Eqs. 117L 1581 . 
The temperature T^, ~ 39 K, see Fig. 8. 



i.e., with the appropriate temperature for the Ohmic 
model, T^^™"^ (for the parameters of Fig. 9, Rh ~ ^o"^ 
at TOhmic ~ 24 K). The interesting feature of Fig. 9 is 

that the dependences H^\t) obtained numerically from 
Eqs. (I2U - GHJ have different behavior near 

Td for T, = 37 K and T,- = 40 K. When T, = 40 K, 
one has Lpb ~ Q everywhere above Td ~ Tj, see Fig. 8. 
In other words, the inductance of the edge wires can be 
disregarded completely. In this case hP{T) has a non- 
monotonic behavior near Td. When Tj = 37 K, one has 



Td w 38.5 K, and this Td lies in the temperature region 
where the inductance of the edge wires is essential. In 
this case the curves h'^\t) of Fig. 9, in fact, monotoni- 
cally tend to zero in the vicinity of Td- As the boundary 
of the temperature region where the inductance of the 
edge wires begins to play a role, one can take the point 
at which Reze — Imze (i.e., at which tpb ss 7r/4). This 
point is just the temperature Tx introduced for the Ohmic 
model. 

In fact, expressions ((5^ . lf37|) are applicable only for 
T > Tx when « 0, and the dependence (T) is 
determined by Ic{T). In the opposite case, at T < the 
phase shift (pi, rapidly tends to tt/2, and expressions H55|l . 
(I57|l become small. This means that the approximation 
used in deriving these expressions fails in this region of 
temperatures. In other words, one cannot neglect the 
time intervals when changes from Ic to —Ic or vice 
versa. However, numerical analysis shows that near Td, 
even though Td < T^, the dependence Hz^\t) obtained 
from Eqs. (ETJ - (EHJ), ^7^, igHI) can be well approximated 
by formula (|17|l if one takes 

(T)«/,(T), (58) 

see Fig. 9. This statement remains true if instead of 
Eq. (|50|) one takes some other form for the temperature 
dependence of Ic- It is only important that the function 
Ic{T) be sufficiently sharp near the point Tj. This find- 
ing means that investigation of Ili^\T) near Td can give 
information on the temperature dependence of the crit- 
ical current Ic in this temperature region, see Appendix 
B. 



VI. CONCLUSIONS 

To describe the magnetic fields generated by the ac 
current in the strip, we formally consider the strip as a 
set of the two edge wires and of wires that form the bulk 
of the strip. The resistances of all these wires can be 
nonlinear functions of the currents flowing in them, and 
the resistances of the edge wires Ri , Rr differ from the 
resistances of the bulk ones. The resistances of the edge 
wires characterize the edge barrier, and they are generally 
different {Ri ^ Rr). To describe the electrodynamics of 
this strip, we also take into account the inductances of the 
wires. Equations (|21|) - 126|l provide the description of the 
magnetic fields and currents in the strip in this general 
case. We then consider in details three special models 
which can be useful for understanding the experimental 
data in various superconductors. The obtained results 
can be summarized as follows: 

Linear ( Ohmic ) model. In this case the resistances of 
all the wires are assumed to be independent of the cur- 
rents flowing in them, and hence the electric fields in the 
strip are linear functions of these currents. To describe 
the known experimental data, it is necessary to assume 
that at Tc the resistance of both edge wires Re is larger 
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than the bulk resistance Rb, but with decreasing tem- 
perature T, Re decreases sharper than Ri,. The temper- 
atures Tsb and defined in Fig. 1 can then be found from 
the conditions: Rb{Tsb) = Re{Tsb), and Rb(Td) « ^J^oUJ. 
The latter condition means that in the Ohmic model at 
Td the ac magnetic field is expelled from the strip due to 
the skin effect, and in the vicinity of this Td the induc- 
tance of the bulk wires becomes essential. Since Re{T) 
decreases sharper than Rb{T) with decreasing T, the in- 
ductance of the edge wires begins to play a role at a tem- 
perature Tx that is higher than Td- This is defined by 
Ileze{Tx) ~ ImZe(Ta;) where Ze is given by Eq. (|35() . 

We now recite the main features of the Ohmic model. 
The first harmonic of the ac magnetic field is practically 
independent of the asymmetry of the edge barrier, rj, 
which is defined by the asymmetry of the resistances Ri 
and Rr, Eq. H37|l . while the second harmonic depends 
essentially on this asymmetry and exists only at 7y 7^ 0. 
The amplitude iji^"* of the first harmonic measured in 
phase with the applied ac current depends on the fre- 
quency Lu only near Td, while the magnitude of the sec- 
ond harmonic of the ac field changes with to only near 
Tx ■ At Tx the magnitude of the second harmonic sharply 
decreases, and a first harmonic that is out-of-phase with 
the ac current appears. A characteristic feature of the 
Ohmic model is also that the ac magnetic field measured 
in units of lao/w, and hence the temperatures T^b, T^, 
Td, are independent of the amplitude lao of the applied 
current. 

Model with nonlinear R^. In this model Re (edge) is 
assumed to be a nonlinear function of the current le flow- 
ing in the edge wires. Namely, we assume the following: 
If le exceeds a critical current characterizing the edge 
barrier, the edge wires do not differ from those in the 
bulk. But this critical current is implied to increase with 
decreasing T, and when it exceeds le, the resistance Re 
sharply drops. Thus, in contrast to the Ohmic case, in 
this model a sharper decrease of Re{T) than of Rb{T) 
is explained by the nonlinear dependence of Re on le- 
The temperature T^b is now near the temperature at 
which the above-mentioned critical current reaches le- 
The temperature dependence of the first harmonic H^^"^ 
is qualitatively similar to that of the Ohmic model, and 
it does not change essentially with the asymmetry of the 
edge barrier. On the other hand, the second harmonic of 
the ac magnetic field is highly sensitive to the asymme- 
try as in the Ohmic case, but now the temperature where 
this harmonic sharply drops is due to the temperature Ti 
at which some characteristic current of the edge barrier 
reaches lao- In this nonlinear model Ti is independent of 
the frequency iv, and it plays the role of the temperature 
Tx marked in Fig. 1. 

Model with nonlinear Rb- In this model Rb (bulk) is 
assumed to be a nonlinear function of the currents flow- 
ing in the bulk of the strip, namely: When the critical 
current for bulk pinning, Ic{T), which increases with de- 
creasing T, exceeds the bulk current /;,, the resistance Rb 
sharply drops. Within this model the temperature Td is 



not due to the skin effect as it occurs for the Ohmic case, 
but now Td is close to the temperature Tj at which the 
critical current in the strip Ic{T) reaches the amplitude 
laO, Ic{Tj) = laO- Thus, the characteristic feature of this 
model is that Td depends on the applied amplitude lao 
of the applied current rather than on its frequency lu. 
The temperature dependence of the first harmonic 
near Td is mainly determined by the dependence Ic{T); 
this enables one to extract information on the critical 
current in this temperature region from the measured 
H^\t), see Fig. 9 and Appendix B. Interestingly, the 
shape of Hz^\t) in the vicinity of Td is different for 
the cases Td > T^ and Td < T^ where T^ is defined by 
KezeiTx) « liRZeiTx), i.e., it marks the point where the 
inductance of the edge wires becomes important. 

Different superconductors generally have different de- 
pendences Rb{T), Re{T), Ic[T), and thus are character- 
ized by different values of the parameters that are essen- 
tial in the formulas of this paper. These parameters also 
depend on the dimensions of the sample and frequency u). 
Therefore, various situations may be expected to occur 
in experiments; compare, e.g., the data of Refs. Q and 
U Possibly, some of them will not be well described by 
the simple models considered above. However, the above 
results allow one to understand the data for a specific 
experiment first qualitatively and then to describe them 
quantitatively by an improved model. 
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APPENDIX A: ASYMMETRIC NONLINEAR 
MODEL 

Here we present formulas for the first and the sec- 
ond Hz harmonics of the ac magnetic field in the frame- 
work of the nonlinear model of Sec. IV in the asymmetric 
case when Ilnl^tx ^ Rex /Ren > 1- Let us introduce the 
following notations: q = w/r, 

I^,^Tex-{l+^^+^y Ia,^IL-{2 + q), 

and 

Ia2^I^ex+IL ■ (1 ' + ' + 

ff Ro/Ren > I en /I ex: wllilc if Rf) / Ren < I en /I ex: WC 

define Ia2 and Jas as follows: 

T TC . . J TC (t , Ren , Ren\ 

Ia2 = lex 1+9 + 9^ , Ia3= I^n 1 + ^ +9-^ ' 

\ i^enj \ iib J 

We also introduce = min(arcsin[/ai//£io], 7''/2) with i =^ 
1-4. 
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The first harmonic Hi^\x) is determined by Ie^\ see 
Eq. (|17|l . This le^^ is described by the following long but 
straightforward expression: 



2/aO Rb 



TT Rt + Re 



sin 2(f>] 



4/c 

— 2£(COS01 



CTi COS 03 - (72 cos 02 j 



-((Tl COS( 



COS 04) + 

sin 201 

2/aO 



2/aO -Rfc 



)2 +0-2 COS 03 

sin 202 



TT Rb + Re 



Rb 



(cos 01 — COS 02) 



TT Rb + Ren 

sin 203 sin 204 
-03 H — 



7r(l + (z) 

4/L . , , , 4/aO /TT 
-(cos 03 —COS 04) H rl — 



7r(l + q) 
sin 204 



2/, 



(72- 



aO 



7r(2 + 9) V2 

+ Ren 



2 y TT q{Rb + Ren) + Re 

sin 202 sin 203 ~ 



(Al) 



where (Ti = 1 if Ro/Ren > I en ex (71 = otherwise, 
while (72 = 1 — (7i. In the symmetric case when Ign/^ex — 

Rex /Ren = 1, OnC haS Ren = Rex — 2i?e, I en = ^ex = 

/^/2, (72 = 0, lai = /a2 = h, = /a4 = ^2 whcrc the 
currents /i, /2 are defined in Sec. IV. In this case formula 
reduces to Eq. (|^ . 
According to formula (|18|l , the second harmonic is de- 



termined by A/i^\ This /S.lY'' is given by 



(2) 



1 cos 01 COS 301 



{Rex + Ren) Rb+Re \3 2 ' 6 

Rb /cos 01 cos 301 cos 02 



Rb+Re 



■(■ 



COS 302 



6 2 
2i?b + i?e„ / sin 202 sin 20i 

' Rb + Ren V 2 2~ 
■OA ■ OA \ , ^'^0 /COS 04 

sm203 -sm204) +JYTqj\~2~ 

cos 304 cos 03 , COS303^ , „^Jln-Jtx),-oM 

- (7 1 [SIW Z02 



(2 + (z) . 
2(1 + (7)^ 



6 / " ^ 2 

qRb - Ren /COS 02 



6 2 

sin 2 03) +(72 /aO /n , 7, N , p V O 
q[Rb + Ren) + Ren ^ 2 

COS 302 COS 03 COS 303 ^ 



6 



6 



(A2) 



(2) 

In the symmetric case the quantity Ale vanishes, as it 
should be. 

We now present formulas for the resistance R of the 
strip in this model. Such formulas may be useful to com- 
pare data obtained from the ac experiments and from the 
conventional transport measurements. Let us define four 
temperatures Ti as temperatures at which Iai{Ti) — 1, 
where Iai{T) = IaiiT)/Iaa, i = 1 - 4, and lao is the 
applied transport current. Then, one has 

RbRe 



R 



Rb + Re 



for T < Tl, 



R = J^'^T (1 - 4) forri<T<r2, 



Rb + Re 
R = aiRb{l-Iex-Ien) 



a2qRbRe 



q{Rb + Ren)+R en 

for T2 < T < T3, (A3) 



R — Rh 



R — Rl 



(1 + 9) 

q 

(2 + g) 



(1 - /,^J for T3 < T < Ti, 



for Ta < T, 



where /g„ = I^^/Iao and /^^ = Note that R is 

independent of the applied current /qq at T < Ti and 
at r > T4. Besides this, if Ro/Ren — Itn/^x some 
temperature , and if the interval T2 — T3 lies above T^ , 
a current-independent resistance occurs in this interval 
as well. 



APPENDIX B: EXTRACTION OF I^{T) FROM 
EXPERIMENTAL DATA 



If the nonlinear Rb model is applicable to describe 
some experimental data on H^\t) near the tempera- 
ture Td, the critical current Ic(Td) can be approximately 
estimated from the relations Ic(Tj) = lao and Td w Tj. 
Changing lao, one finds Ie{T) in some temperature inter- 
val which is determined by the range of lao ■ On the other 
hand, the temperature dependence Ic{T) can be also es- 
timated from the shape of the curves H^\t) measured 
at fixed lao- Using the approximation H58|) and formula 
(I17|l . we obtain 



Ic{T) 



2wHPiT)-Iao-g 

f-9 



(Bl) 



where H^\t) is the first harmonic measured at an in- 
ner point X (i.e., at \x\ < w), g = {xw / it) / (iiP' — x^). 



/ = 2wHP{Te)/Iao, and H^z''{Te) is the first harmonic 
H^^^ measured at the same point a; at a temperature near 
Tc- Here, to express / in terms of H^^'' (Tc), we have used 
formula (jl7|l again and have taken into account that at 

T = Te the edge barrier is absent, and I^^^ = lao- When 
the sample is an "ideal" strip (i.e., uniform with constant 
thickness, width, and homogeneous critical current den- 
sity), one has / = (l/27r) ln(|w — x\/\w + x\). Note that 
these two ways of estimating le more reliably give the 
temperature dependence of le rather than the absolute 
value of the critical current. Indeed, a,t T — Td when one 
has Hi^\T) — 0, the first method gives Ic ~ laO, while 
the second method leads to Ic — laO ■ 9/ (.9 — /) ~ ^ao/2. 
This discrepancy is caused by the inaccuracy of the rela- 
tion Td « Tj and of the approximation (|58|l . 
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